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Abstract. Relativistic electrons, seen in the large diffuse 
radio halo of the Coma cluster of galaxies, should scatter 
background photons to higher energies. We calculate the 
inverse Conipton contributions from the microwave back- 
ground, from the local radiation field of elliptical galaxies 
and from the thermal X-ray emission of the intra cluster 
medium to the different observed energy bands, and draw 
restrictions to the shape of the spectrum of the relativis- 
tic electron population. The expected electron spectra in 
different halo formation models are discussed, and signa- 
tures for a future distinction of these models are presented, 
tracing the injection processes and the influence of optical 
thickness of the radio halo for very low frequency emission 
due to synchrotron self absorption. Some of the upper lim- 
its on the electron spectrum translate into lower limits of 
the halo magnetic fields, since the observed synchrotron 
radiation is known. Despite uncertainties in the extrapo- 
lation of the radio spectrum to lower frequencies a lower 
limit to the central field strength of Bo > 0.3 seems 
to be sure. This is comparable to the value resulting from 
minimal energy arguments. If Bo < 1.2 /iG, the recently 
reported extreme ultraviolet excess of Coma (Lieu et al. 
1996) could result from inverse Compton scattered mi- 
crowave photons. Formulae for transrelativistic Thomson 
scattering are given in the Appendix. 

Key words: magnetic fields - galaxies: clusters: individ- 
ual: Coma - (galaxies:) intergalactic medium - radiation 
mechanism: non-thermal - radio continuum: general 



1. Introduction 

An estimate of a lower limit on the magnetic field strength 
in the Coma halo region by Rephaeli et al. (1994) is 
B > 0.1 /iG, using the radio spectra of the halo and 
an upper limit to inverse Compton (IG) photons in the 



hard X-ray region measured by the OSSE experiment. The 
limit on the field strength derived from the upper limit to 
the gamma-ray flux above 100 MeV by EGRET resulting 
from relativistic bremsstrahlung is i? > 0.4 fiG (Sreeku- 
niar et al. 1996), if the radio spectrum can be extrapolated 
to lower frequencies. The field strength, derived by Kim 
et al. (1986), resulting from minimal energy arguments 
is i? « 0.6/iG(l -I- fcp)^/^/i5Q^, where fcp is the proton 
to electron energy ratio. Measurements of Faraday rota- 
tion of polarized radiation seen through the Coma cluster 
medium, combined with gas profiles, derived from X-ray 
observations of the hot intra cluster medium (ICM), give 
magnetic fields of 1.7 fiG (Zfts/lO kpc)^^/^ /i^q^ (Kim et al. 
1990) and 6.0^*0 {lhs/l]<ipc)-^^^hl()^ (Feretti et al. 1995). 
Ifis denotes the field reversal scale of the magnetic field 
and can in principle be measured by depolarization obser- 
vations. In both measurements the resolution of /frs was 
limited by the resolution of the telescope, and therefore 
higher field values could result from a smaller field rever- 
sal scale. Strong magnetic fields on the order of a few tens 
of /xG are expected from the injection of radio plasma from 
radio galaxies into the IGM (EnBlin et al. 1997). Although 
field limits derived by IG flux limits are lower than the 
rotation measurements, they have the advantage of being 
less independent of the model. Quantities in this article are 
scaled to a Hubble constant of Ho = 50 km s~^ Mpc~^ /i, 
where h^o indicates their scaling. 
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2. Multifrequency Observations 

2.1. Radio 

The radio flux spectrum of the diffuse radio halo Coma G 
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(Kim et. al 1990) results from a relativistic electron pop- 
ulation, which we assume to be 



N,{r, /5e7o) d(/3c7c) = Ce(r) (/3e7c)"" d(/3c7c) 



(2) 



with a spectral index a — 3.68±0.12. Pc — /3c 7o Wc c is the 
electron momentum, and /3c should not be confused with 
the cluster shape parameter (3. The synchrotron flux of a 
power law electron distribution in an isotropic distribu- 
tion of magnetic fields within the halo volume (Eq. (6.36) 
in Rybicki & Lightman 1979), averaged over an isotropic 
distribution of electron pitch angles, is 
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Bo, Ceo are taken at the cluster center, D = 139Mpcft,gQ 
is the distance to Coma, F is the gamma function, and 
Vo — 3eBo/{2TTmcc). We assume that magnetic field en- 
ergy density and relativistic electron density scale radi- 
ally with the same /3-profile as the background gas does, 
- (1 + (r/rc)2)-3/3/2 with 13 = 0.75 and = 400 kpc 
(Briel et al. 1992). This is reasonable since every possi- 
ble source of magnetic fields such as injection by radio 
galaxies, compression of primordial fields frozen into the 
plasma, and amplification of fields by turbulent gas mo- 
tion should result in a rough scaling of the magnetic en- 
ergy density with the thermal energy density. Similar ar- 
guments hold for the supply of the relativistic electron 
population. But the best support of this assumption is 
given by Fig. 3 of Deiss et al. (1997), which compares the 
radial X-ray profile with the radio profile of the cluster, 
and shows that the latter is steeper. The X-ray emission 
scales with the gas density as ~ (r) , and the radio emis- 
sion as ~ Cc(r) i?^"+^^/^(r) ~ ng-^^(r), and is therefore 
a little bit steeper, if we use the scaling assumed above 
(i?^(r) ~ Cc{r) ^ ne{r)). The resulting projected emis- 
sion profile is similar to that measured by Deiss et al. 
(1997). The emission profile itself is shown in Fig. |^. 

The emission weighted volume within the halo radius, 
which is at least R = 1.2 Mpc /i^q^ (Deiss et al. 1997), is 
given by 



Vsync = 2TT^BrlB ^2 (3,3 (S±5^_l)) ^ 



(5) 



where Bx{a, b) denotes the unnormalized incomplete beta 
function (Eq. (6.6.1) Abramowitz & Stegun 1965), and 
should not be confused with the symbols for the field 
strength B{r) and Bq. (j>B is the filling factor of the mag- 
netic field in the volume occupied by the relativistic elec- 
trons, which we expect to be close to unity. Comparing 
Eq. |l| and Eq. 3 gives 



Ceo = 3.28 • 10"^ cm"3 h^Q 



Bo 



-2.34 



2.2. Infrared 

Wise et al. (1993) searched extensively for diffuse far- 
infrared emission from clusters, and detected Coma 
marginally at 60 /im, but not at 100 /xm. Their estimate of 
the excess fluxes above the background within 30 arcmin 
of the cluster center are 5 ± 4 and 1 ± 10 mJy arcmin"^. 
(1 arcmin = 39 h'^Q kpc). We use two sigma upper limits 
in the following, namely 13 and 21 mJy arcmin"^. 

2.3. Extreme Ultraviolet 

Recently, Lieu et al. (1996) have reported the detection 
of extended, extreme ultraviolet emission from the clus- 
ter center with the Extreme Ultraviolet Explorer (EUVE). 
They interpret this as evidence for relatively cold gas com- 
ponents with temperatures of kT — 0.07 and 0.4 keV. But 
Dixon ct al. (1996) tried to detect resonance line emis- 
sion with the Hopkins Ultraviolet Telescope (HUT), which 
would be expected in this case. Their nondetection does 
not reject the thermal gas model completely, but gives 
additional constraint on the properties of the cooler ICM 
components. The EUV emission might also arise from in- 
verse Compton scattering of microwave or starlight pho- 
tons, giving then direct informations about the relativistic 
electron population. On the other hand, if the real nature 
of this emission can be proven to be different from IC, 
then we can interpret their flux as an upper limit to any 
possible IC contribution, and therefore giving an upper 
limit to the electrons. If the detailed physics of such an 
emission process can be understood, this limit could be 
improved. 

Their EUVE observation covered a field with a radius 
of 30 arcmin centered on the X-ray center of Coma, and 
they detected emission within the central 15 arcmin. The 
count rate in the passband of 65 to 245 eV was 36 % above 
the expected count rate from the emission of the thermal 
gas with kT = 8.2 keV. We estimate the latter flux in 
this band to be 8.1 • 10~^cm~^s~^, using their emission 
measure, which is in good agreement with that expected 
from the /3-model of the gas with a central gas density of 3- 



10 



'3cm~3 h^Q^^, and using the Bremsstrahlung formulae 
given in Novikov & Thorne (1972). This translates to an 
excess flux of 2.9-10^^ cm^^ s~^, which we use as an upper 
limit to the IC flux in the 65 to 245 eV passband in the 
following. 

2.4. High Energy X-Rays and Gamma-Rays 

Rephaeli et al. (1994) measured the high energy X-ray 
(HEX) flux from Coma with OSSE to be lower than 
6 • lO^^cm-^s-ikeV-i within the 40-80 keV band. The 
EGRET flux limit on gamma rays from Coma is [E^ > 
100 MeV) < 4 • 10"^cm-2 s'^ (Sreekumar et al. 1996). 



(6) 3. Inverse Compton Emission 
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Bockground Photon Fields: 

MWB= microwave background 
opt = optical emission 

from elliptical galaxies 
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Fig. 1. The central electron population ol Coma. The line labeled with 'Bremsstrahlung' shows the thermal population. Also 
limits from bremsstrahlung in the OSSE and EGRET bands are shown. The other lines result from radio observations and 
are labeled with the corresponding central field strength. Limits from inverse Compton scattering are labeled with the relevant 
photon population and the resulting energy range. All limits are two sigma limits except the EUV limits, which result from a 
detected flux, and might be read as data points, if the EUV emission results from inverse Compton scattering instead from a 
cool ICM component. A horizontal line in this diagram corresponds to a power law Nc{Pc) ^ Pc ■ 
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Table 1. The IC limits on the normalization constant Ceo of the power law electron distribution (Eq. g) resulting from the 
different combinations of background photon field and observation bands, and the typical electron momenta necessary for this 
scattering. 



3.1. Coma's Electron Population 

Only two regions of the electron spectrum are visible: 
the X-ray emitting thermal bulk with a temperature of 
8.2 keV and relativistic electrons around 1 GeV visi- 
ble in the radio. Outside of these regions information 
can be obtained from the detection or nondetection of 
IC flux resulting from electrons and background photon 



fields, mainly the microwave background (MWB) and the 
starlight of the cluster members. The radio photons of 
the halo do not contribute significantly to the IC flux: 
The energy density of the radio emission is of the order 
of iradio/(47r cr^) « 10~^ eV cm~'^ and therefore 6 orders 
of magnitude lower than that of the MWB. The limits re- 
sulting from the different combination of photon fields and 
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observation bands are given in Tab. |l| and plotted in Fig. 
|l| above the typical electron momenta necessary to scatter 
the peaks of the blackbody spectra into the observation 
bands from the relation < e >= | 7^ 2.70 kT (Blumenthal 
& Gould 1970). For the calculations, which are explained 
below, the spectral index of the radio electrons is used. 

3.2. Microwave Background 

The expected IC flux from the MWB is (derived from 
Eq. (7.31) in Rybicki & Lightman (1979) for Thomson 
scattering) 



/ic(a) = 



A,^„r,V,c«.)c„(§| ,7, 



+ Aa + ll) 



(a + 3)2(a + 5)(a2 
Vic = 2nrlB ,2 (|,|(/3~1)) , 



^r(^)C(^) (8) 

(9) 



where C denotes the Riemann zeta function, and Rq the 
radius covered by the observation. The electrons seen in 
the radio scatter microwave photons into the hard X-ray 
band. Thus, with the help of Eq. ^ we can use their limit 
in order to get a lower limit on the central field strength of 
Bo > 0.2 fiG (j)g'^'*^h^Q, which is independent of the Hub- 
ble constant. Assuming a uniform magnetic field strength 
over the whole cluster volume by setting /? = 0, gives only 
Bo > 0.1 fJ,G 4)g'^'*^hc,Q, in agreement with the value of 
Rephaeli et al. (1994). But this is an unrealistic configu- 



ration, as explained in Sect. 2.1 



Since none of the halo formation theories discussed in 
Sect. ^ predicts any break in the electron spectrum be- 
tween 100 MeV/c and GeV/c, it is reasonable to extrap- 
olate the spectrum and use the limit given by the EUV 
flux to predict a field strength stronger than 



Bo > 1.2 ^G 



i>B FEUV.IC-limit 

0.16/(cm2skeV) 



-0.43 



(10) 



If the EUV IC contribution can be further constrained, 
a higher field strength would follow. A break by 0.5 in 
the radio spectral index, placed at the lowest observed 
frequency of 30.9 MHz (Henning 1989) still gives a limit 
Bo>0.i^iG, 



-0.54 



3. 3. Optical Photons 

The luminosity of elliptical galaxies within the central 
700arcmin2 is given by an integration of the R-band lu- 
minosity function of Seeker & Harris (1996) plus the lu- 
minosities of NGC 4874 and NGC 4889 (Strom & Strom 
1978), which are not included in this luminosity function. 
We assume that this radiation has a blackbody spectrum 
with a typical temperature of 3000 K, and therefore use a 
bolometric correction of tor — mboi — 1.3 (Webbink & Jef- 
fers 1969). The radial emission profile is that of the galaxy 



distribution. e{r') - (l + (r7rG)^)""^', with = 0.8 and 
rQ — 160kpc/i^o^ (Girardi et al. 1995), which we use up 
to a radius of Rq = 5Mpc/i^Q^. We estimate a central 
bolometric emissivity of So — 1-0 • lO^'^ Mpc~^ /150. In 
order to get the photon density field at a position r, this 
has to be folded with the radiation field of a source at 
r': n^{r,r') — e(r')/(47rc|r — r'p). Integrating over all 
possible angles between r and r' gives the photon density 



Rg 



r -\- r' 



(11) 



Normalizing this by dividing with the photon density of a 
blackbody cavity, and integrating it together with the nor- 
malized /3-profile of the relativistic electrons numerically 
over the volume gives the effective volume Vic, which re- 
places Eq. 9. The scattering of optical photons into the 
EGRET band needs electrons, visible in the radio. Thus 



a limit of Bo > O.lfiGc 



-0.43 



can be set. The assumed 



temperature of the radiation mainly enters this calcula- 
tion by the bolometric correction. The anisotropy of the 
optical photon field at a given position does not change 
the resulting IC flux of the cluster, since the integration 
over shells of constant radius averages the different con- 
tributions. 

3.4. Thermal and Nonthermal X-Ray Photons 

The radial density profile of X-ray photons is also de- 
scribed by Eq. with an emissivity profile e(r') ~ 
n^{r') ~ (1 -I- {r' /rc))~^^ . We calculated the spectrum 
of the Bremsstrahlung using the Gaunt factor given in 
Novikov & Thorne (1972), a central gas density of n^o = 
3 • 10"'^ cm~'^ hgQ^, and a gas temperature of kT = 8.2 
keV (Brief et al. 1992). The IC scattering of this photons 
had to be calculated numerically. The contribution to the 
EGRET band is negligible, due to the reduction of the IC 
cross section in the Klein-Nishina regime. The calculation 
of the scattering into the OSSE band was done using the 
exact probability distribution of frequency shifts of Thom- 
son scattering, valid from non- to ultrarelativistic electron 
momenta, which is derived in the Appendix. 

The limit resulting from the expected bremsstrahlung 
in the OSSE band of a suprathermal electron distribution 



is Ceo < 6 ■ cm^-^ h 



3^,1/2 



50 



Bremsstrahlung of relativistic 
electrons was calculated with the formulae in Blumenthal 
& Gould (1970). The EGRET limit translates into Ceo < 
3.3-10-2 cm /igg . These limits are also shown in Fig ^ 



4. Radio Halo Formation 

The slope of the electron spectrum is determined by 
the acceleration, cooling and injection processes. Three 
sources for the energetic electrons have been discussed in 
the literature: 
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— The primary electron model (Jaffe 1977, Rephaeli 
1979), in which the electrons diffuse out of radio galax- 
ies and form the radio halo. Their spectrum should 
be a straight power law. A steepening of the spectral 
index appears at a break frequency, if there is some 
momentum dependent escape from the halo region. 

— In the in-situ acceleration model (Jaffe 1977, Roland 
et al. 1981, Schhckeiser et al. 1987, Burns et al. 1994, 
Deiss et al. 1997) the electrons are accelerated within 
the ICM by plasma waves or shocks. Burns et al. (1994) 
report evidence for a recent merger event of Coma ^ 2 
Gyr ago, which was strong enough to power the ra- 
dio halo. An exponential cutoff in the spectrum should 
occur at the momentum, where the cooling time scale 
gets comparable to the acceleration time scale. 

— The secondary electron model (Dennison 1980) as- 
sumes the relativistic electrons, but also positrons, to 
result from pion decay after hadronic interaction of 
a relativistic proton population with the background 
gas. Since the protons should result from radio galax- 
ies (Enfilin et al. 1997) and supernova remnants in 
starburst galaxies (Volk et al. 1996), which are able 
to accelerate to high energies, no cutoff in the elec- 
tron spectrum is expected. But the injection spectrum 
peaks at 35 MeV/c (due to the kinematics of the decay 
TT^ v^jjL^^lx^ Ufj^VeC^), and therefore the equilib- 
rium spectrum should flatten to a spectral index a — 2 
(due to synchrotron and IC cooling) below the main 
injection energy range. 

Schhckeiser et al. (1987) report a high frequency cutoff 
in the radio spectrum of the halo at 2.7 GHz, which would 
indicate a cutoff in the electron distribution, and therefore 
supports the in-situ acceleration model. But Deiss et al. 
(1997) claim that their much higher measured flux at 1.4 
GHz is not consistent with the value of Schhckeiser et al. 
(1987). Deiss et al. (1997) suspect that the value of the 
flux given by these authors is much too low, which im- 
plies that the claimed strong steepening of the spectrum 
above 1.4 GHz is not real. In the literature (e.g. Schhck- 
eiser et al. (1987) or Deiss et al. (1997)) citations of an 
upper flux limit at 5 GHz appear frequently, measured by 
Waldthausen (1980). But his observation was restricted to 
the central region of the Coma cluster, and therefore no 
proper background determination was possible. In order 
to clarify the question about the existence of a cutoff in 
the electron distribution, which is important for models 
of halo formation, new high frequency observations of the 
diffuse radio halo are necessary. 

The low energy spectral break expected in the sec- 
ondary electron model cannot be used easily for a distinc- 
tion between the models, since we expect a sharp step in 
the spectrum in the same energy range, which could make 
it difficult to detect the appearance of the break. This step 
is caused by a sharp change in the cooling time scale as 
a function of energy, due to the transition from optical 
thin synchrotron emission at higher energies to optically 
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T(r)/(«.(0) 
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Fig. 2. Normalized optical depth r(r) for synchrotron self ab- 
sorption (solid) and synchrotron emissivity Svir) (dotted) as 
a function of radius. The numbers indicate the percentage of 
total emission produced within the indicated radius. A radial 
unlimited emission region was assumed for simphcity. 

thick emission below the step energy. In the optical thick 
regime, the radio power emitted by the electrons is reab- 
sorbed, and therefore synchrotron cooling is suppressed. 
The synchrotron self absorption coefficient (Eq. (6.53) in 
Rybicki & Lightman 1979) averaged over an isotropic dis- 
tribution of electron pitch angles is 
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and can be evaluated at the cluster center using Eq. ^: 
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(14) 
(15) 



where the frequency was translated to electron momentum 
using the monochromatic approximation v = (/3c 7c)^ i^o- 
The optical depth T(r) at a radius r is given by an integral 
from that radius to infinity over the radial profile of the 
absorption coefficient 



C,{r)B{r)^ 



[l + {r/r,r 



i/3 



(16) 



and is plotted in Fig. ||. Most of the total power at a given 
frequency is emitted within 2 core radii, and would result 
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from a region with optical depth bigger than one if the 
central absorption coefficient is ai,(0) > 1/(0.05 rc). This 
happens for i' < D — 21 kHz (i3o/MG)°'^'^ correspond- 
ing to Pe < f'c = 81MeV/c (Po/mG)"°-^^ This sudden 
change in the radio spectrum can hardly be detected, but 
should make a signature in a halo electron population, 
which was injected at higher energies and cools by IC and 
synchrotron losses, as it is the case for the primary and 
secondary electron model. The time independent solution 
of the continuity equation for this electrons in momentum 
space is 



iVe(Pe) 



1 



dpqdp) , 



(17) 



where qdp) is the injection spectrum and 



6(Pe) = |^T f-^ 

o yrUc c 



1^ H(Pe - Pe) + C/ph ) (18) 
6 TT 



the cooling function. J7ph is the photon energy density, 
which is dominated by the MWB, and H is the Heavi- 
side step function, used as an approximation of the sharp 
change in opacity. The height of the step is therefore 



(19) 



The exact shape of the transition is a little bit more com- 
plicated to calculate, since the problem is nonlinear due 
to the mutual dependence of the electron spectrum and 
the absorption coefficient, the geometry of the radiation 
transport, which has to be considered, and the time depen- 
dence before establishing the equilibrium spectrum. For a 
field of e.g. Bo ~ 10 /iG a step in the electron spectrum 
of one order of magnitude is expected at Pe = 26MeV/c 
(using our simple approach), very close to the position of 
the expected flattening below the injection peak in the 
secondary electron model. The in-situ acceleration model 
predicts no step at this energy, since the cooling func- 
tion determines only the high energy cutoff region of the 
spectrum, where the cooling gets comparable to the accel- 
eration. We note, that any break in the electron spectrum 
below the observed range would lower the effect of syn- 
chrotron self absorption. Assuming a break in the radio 



spectrum at 30.9 MHz as in Sect. 3.2 shifts the step in the 
above example to 15 MeV/c. 

5. Conclusion 

The search for diffuse emission in any photon energy band 
leads to interesting information about the shape of the rel- 
ativistic electron population in the Coma cluster of galax- 
ies. We present a compilation of several limits, and discuss 
spectral features of the electrons within the various ra- 
dio halo formation models. A signature in the low energy 



range is expected to arise from the sharp change in opac- 
ity for low frequency radio photons due to synchrotron self 
absorption, which could distinguish between a relativistic 
electron population which is cooling and one which is gain- 
ing energy, as it is predicted for primary and secondary 
electron models on the one hand, and for in-situ accelera- 
tion models on the other hand. Another signature of the 
secondary electron model is a flattening of the spectrum 
below the main electron injection energy. Future observa- 
tions of diffuse emission might test the electron distribu- 
tion sensitively enough in order to see these signatures. 
Present observations restrict the electron population suf- 
ficiently, so that a central magnetic field strength higher 
than 1.2 /iG is needed in order to explain the observed 
synchrotron emission, if the steep volume integrated ra- 
dio spectrum with spectral index of 1.34 can be extrap- 
olated to lower frequencies, or a limit of Bo > 0.3 /xG, 
if a break appears below the observed range. Since these 
limits arise from a detected EUV flirx excess, which is 
interpreted to result from a cool gas component of the 
intra cluster medium, the possibility remains, that this 
flux results in fact from inverse Compton scattering of mi- 
crowave photons, if the field strength is below 1.2 /iG. This 
possibility seems to be in contradiction with the Faraday 
measurement of Feretti et al. (1995) of 6/iG. But since 
the inverse Compton argumentation gives volume aver- 
aged field strength which are weighted with the relativis- 
tic electron distribution, stronger field strength in regions 
with a low density of relativistic electron could solve this 
contradiction. If the nature of this EUV excess can be un- 
derstood theoretically, the EUV limit on inverse Compton 
flux might be improved and a higher field strength limits 
might result. 

Limits without extrapolation of the electron spectrum 
come from the scattering of starlight photons to EGRET 
energies {Bo > 0.1 /iC) and microwave photons to OSSE 
energies (Pq > 0.2 /iC). Magnetic field limits from in- 
verse Compton limits are below the estimates resulting 
from Faraday rotation measurements, being therefore in 
good agreement. They measure the magnetic content of 
the Coma cluster on a more global scale, and demonstrate 
independently that the observed Faraday rotation results 
from magnetic fields within the intra cluster medium. Fu- 
ture EUV, X- and gamma-ray telescopes will raise these 
limits on the magnetic field strength in the Coma halo 
region to the level of Faraday measurements, or they will 
detect energetic nonthermal photons from Coma, resulting 
from inverse Compton scattering or from other sources. 
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A. Transrelativistic Thomson Scattering 

We follow the derivation of the formulae of transrelativis- 
tic Thomson scattering by Wright (1979). Scattering of 
a photon from an isotropic photon field by an electron 
with velocity /3c c and energy 7c rrtc has a probability 
distribution of the angle 9 between photon and electron 
direction in the electron's rest frame given by 



(Al) 



where fj, = cos 6. The logarithmic frequency shift from ly 
to p' of the photon after scattering into an angle 9' is 



In 



In 



1 + /3c// 



(A2) 



1 - /3o/i . 

The probability distribution of fi' — cos 9' for a given fi is 

g{p'; /i) d^l' = ^[1 + AiV" + ^(1 - m')(1 - A*")] dp' .(A3) 

The probability for a shift s given the electron velocity is 

'dp'' 



P{s;Pc)ds 



dpf{p)g{p';p) 



ds 



ds 



(A4) 



p' {s, (3ct p) is given by (A2) and the range of integration 
is given by the conditions —l<p<l and —l<p'<l. 
Rephaeli (1995) simplified this integral by assuming that 
the scattering is isotropic in the electron's rest frame, 
which allows to average over the angle between back- 
ground and scattered photon. He compared scattered pho- 
ton spectra using the simplified and exact formulae nu- 
merically, and found no significant difference for electrons 
with energies of 1...15 keV, typical for the thermal compo- 
nent of the ICM. But suprathermal electrons, expected to 
be the connection between the thermal and the halo elec- 
trons, see an anisotropic photo n fi eld in their rest frame. 
Fortunately, the integrand in (A4) is a rational function 
of p and therefore analytically integrable. Integrating the 
two regimes s < and s > separately and introducing 
the dimcnsionless momentum p = /3c 7o leads to: 



3e"(l-t-e'') 



2p 



iS-\s\) 



+ ^^^ [1 - + 3e^(l - e^)(9 + 8p^ + Ap^)] (A5) 
The maximal logarithmic frequency shift is given by 

'i + A 



s — In 



111 



V^TjP -p. 



(A6) 



and therefore is P{s,p) = for \s\ > s. The distribution of 
frequency shifts of scattering with a spectrum of electrons 
ndp) dp 



p/ X , / dpnc{p)P{s\p) 

Pis) ds = 7^- -j—^ ds 

J dpucip) 



(A7) 



has to be folded with the photon spectrum n(j/) dv in order 
to obtain the source spectrum of scattered photons 



q{h'') dv' = (Jt^ cUc 



+ 00 



dsP{s)n{v' e^') 



dv' . (A8) 



These formulae can be applied from non- to ultrarelativis- 
tic electron momentum, they are valid for scattering to 
higher and lower energies compared with the background 
photon energy, but they are restricted to the Thomson 
regime {hv <^ 7cmoC^). Applications are especially IC 
scattering by transrelativistic electrons, such as suprather- 
mal electrons, or electrons of a very hot plasma. Since the 
evaluation of (A5) needs only a little more computation 
time than its approximations, it could be included in any 
numerical code for Thomson scattering. 
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